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I. 


INTRODUCTION 


In recent years, active and passive control of sound and vibration in aeroelastic 
structures have received a great deal of attention due to many potential applications to 
aerospace and other industries. There exists a great deal of research work done in this 
area. Recent advances in the control of sound and vibration can be found in the several 
conference proceedings. In this report we will summarize our research findings 
supported by the NASA grant NAG- 1-1 175. 

The problems of active and passive control of sound and vibration has been 
investigated by many researchers for a number of years. However, few of the articles are 
concerned with the sound and vibration with flow-structure interaction. Experimental 
and numerical studies on the coupling between panel vibration and acoustic radiation due 
to flow excitation have been done by Maestrello and his associates at NASA/Langley 
Research Center. Since the coupled system of nonlinear partial differential equations is 
formidable, an analytical solution to the full problem seems impossible. For this reason, 
we have to simplify the problem to that of the nonlinear panel vibration induced by a 
uniform flow or a boundary-layer flow with a given wall pressure distribution. Based on 
this simplified model, we have been able to study the control and stabilization of the 
nonlinear panel vibration, which have not been treated satisfactorily by other authors. 

The vibration suppression will clearly reduce the sound radiation power from the panel. 
The major research findings will be presented in the next three sections. In Section II we 
shall describe our results on the boundary control of nonlinear panel vibration, with or 
without flow excitation. Section III is concerned with active control of the vibration and 
sound radiation from a nonlinear elastic panel. A detailed description of our work on the 
parametric vibrational control of nonlinear elastic panel will be presented in Section IV. 
This paper will be submitted to the Journal of Acoustic Society of America for 
publication. 


II. BOUNDARY CONTROL OF NONLINEAR PANEL VIBRATION 

This work is concerned with the analytical study of a stabilization scheme to 
suppress the large amplitude panel vibration by means of boundary damping: that is, by 
introducing extraneous dissipative boundary conditions, an otherwise unstable system 
will achieve its global and exponential stability. To be specific, we consider the panel 
vibration to be governed by a nonlinear beam equation, with or without aerodynamic 
forcing, subject to the clamped-free boundary conditions: that is, the left end of the panel 
is clamped and the right end is free to move. Without any control, the compressive force 
and the aerodynamic loading may be time-dependent. To stabilize the system, we apply 
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to the free edge a set of boundary controls, which consists of the combination of a 
vertical force, a bending moment and a tensile force. By means of the energy method and 
some mathematical inequalities, the boundary stabilization of the model problem for a 
vibrating nonlinear elastic panel was treated analytically. In general, the panel is subject 
to a compressive in-plane loading combined with an aerodynamic forcing. Without any 
control, the panel may buckle and would oscillate due to flow induced instability. To 
stabilize the panel, a set of boundary controls was introduced as the combination of a 
bending moment, a vertical point force and a tensile force applied to the free edge. Two 
cases, corresponding to the absence and the presence of an aerodynamic loading, were 
treated separately. 

Without the flow, this is the case of free vibration. Even though the energy of the 
uncontrolled system is conserved, with an initial disturbance, the system may buckle or 
sustain persistent large amplitude oscillation. To render the energy an exponential decay, 
it was found sufficient to apply a tensile force, if necessary, to reduce the net force to a 
subcritical level and, at the same time, to introduce boundary damping. The damping 
mechanism consists of frictional forces and torque, which are linearly proportional to the 
transverse and the angular velocities of the right edge, respectively. Therefore, if the 
compressive force is subcritical, the passive control in the form of boundary damping 
suffices to stabilize the system. In an analogous situation, the result seems to be in 
agreement with the experimental evidence that boundary damping is effective in 
suppressing the panel vibration. Jn the presence of the unsteady flow and a compressive 
force, the panel may buckle or flutter. If the flow velocity is oscillatory and decays 
rapidly, it is possible to stabilize the panel by applying a time-varying tensile force 
together with boundary damping as before. However, for a slowly varying compressive 
force and, at the same time, with a small flow parameter the system can be stabilized as 
in the free-vibration case. For the purpose of illustration, several numerical examples 
were provided. In what follows, a few comments are in order. 

(1) The paper deals with a model problem which may not be realistic in practice. 
However, at least qualitatively, the analytical results support the experimental 
evidence that the passive control by boundary damping is effective in 
suppressing panel vibration. 

(2) As in the stability analysis of an ordinary differential equation, we set the flow 
pressure equal to zero, because one is interested only in the stability of the 
equilibrium solution due to an initial perturbation. In the presence of the 
pressure fluctuation, the system vibration may be suppressed by a distributed 
control. This will be discussed in the next section. 

(3) The theoretical approach based on the energy method and some mathematical 
inequalities, when applicable, is rather powerful. However, it certainly has its 
limitations. In general one needs both the analytical and numerical techniques 
to tackle the difficult the flow-structure interaction problems. 

The results of this research and related work are summarized in our papers [l]-[3]. 
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III. ACTIVE CONTROL OF SOUND AND VIBRATION 

The main purpose of this work is to study, based on a nonlinear panel model, the 
effectiveness of the active control to suppress the panel vibration and sound radiation 
induced by the unsteady pressure fluctuation. The control consists of a distributed force 
applied normally to one side of the wall. For simplicity, the flexible panel is assumed to 
be hinged to the rigid plates at both ends. The coupled equations governing the nonlinear 
panel vibration and acoustic radiation problem were given. For the optimality criterion, a 
time-average energy or objective functional was introduced to measure the performance 
in controlling the vibration and sound radiation. By applying the variational method, we 
derived the optimality equation for the control force distribution which is coupled with 
the controlled equations of motion. By using an eigenfunction expansion, the modal 
control problem was formulated. The truncated modal control problem was solved 
numerically by the shooting method for a two-point boundary value problem in the time 
domain. The numerical results were obtained to demonstrate the effectiveness of active 
control. 

The main results of this research are summarized as follows. 

(1) For the control of panel vibration, given the control objective function, the 
optimal control can be found in the form of an external pressure applied to 
the wall. We derived the partial differential equation for the optimal control 
which is coupled to the equation of motion. 

(2) For the given objective function, it is possible to derive the optimality system 
for the control of sound radiation governed by the wave equation. 

(3) In both cases, the optimality system consists of a coupled nonlinear boundary 
value problem in space and time. 

(4) In the case of truncated modal control, the optimality system yields a two- 
point boundary value problem for a finite set of nonlinear ordinary 
differential equations. 

(5) The truncated optimality system for modal control was solved numerically. 
The results show the following: (i) for linear panel vibration, the control is 
highly effective and can almost completely eliminate the vibration over a 
short time horizon, (ii) In contrast, a non-linear panel, in general, responses 
less sensitively to the active control. The control is more effective at lower 
vibration frequencies and with weak nonlinearity. For fixed nonlinearity, the 
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effectiveness of control diminishes as the frequencies increase, and 
eventually the system loses control completely, (iii) By applying the optimal 
vibration control, the sound radiation intensity can also be reduced 
significantly. 

The research is related to the numerical studies performed by Maestrello and 
some of his collaborators [1], [2], The results of our work can be found in the 
published paper [3], 
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IV. VIBRATIONAL CONTROL OF NONLINEAR PANEL 
1. Introduction. 


The problem under consideration is the stabilization of the nonlinear panel oscillation 
by an active control with a vibrational actuator. This work was motivated by the recent experimental 
investigations of the second author (Maestrello [1], [2]), who demonstrated clearly that the 
vibrational control could be an effective means of stabilizing the boundary-layer flow as well as the 
panel vibration. This paper will offer a general method of vibrational control and its application to 
the problem involving a nonlinear elastic panel excited by the periodic wall-pressure fluctuation in a 
boundary-layer flow. 

The general principle of active control, the vibrational control in particular, is to introduce an 
action which affects a change in the behavior of a dynamical system in a desirable manner. In the 
boundary layer transition control [1], the periodic heating and cooling of the wall induce a 
parametric vibration of the fluid viscosity which, in turn, stabilizes the flow. In the case of panel 
vibration [2], a properly added vibrational force with the same forcing frequency may result in 
suppressing the subharmonic oscillations, (see section 2). The suppression of subharmonics, of 
course, has the implication of controlling the chaotic motion. 

The main idea of vibrational control stems from the fact that an inverted pendulum can be 
stabilized at its upper equilibrium position when the lower suspension point executes a rapid vertical 
vibration, (see e.g. [3]). Based on this idea, a general principle of vibrational control was proposed, 
notably by Meerkov [4], to stabilize the equilibrium points of some finite-dimensional linear 
systems. Application of this principle to reactor dynamics was done by Bellman et al [5], By 
contrast, in this paper, we will extend this control principle to stabilize the periodic motions of 
infinite-dimensional systems, instead of equilibrium points in finite dimensions. In addition to the 
high-frequency parametric vibrational control used in [4,5], a vibrational force with the same forcing 
frequency will be required. Unlike the usual feedback of feed-forward control, the vibrational 
control does not need accurate measurement of the system inputs and outputs and can be 
implemented much more easily, especially for an infinite-dimensional system under consideration. 

In this paper we consider the panel vibration which satisfies the initial boundary value problem 
for the nonlinear beam equation [6]: 


r 

< 


mdjw + cd,w - [Q + N(t)]d 2 x w + Dd x w = Ap(t,x),0 < x < l 
vv(f,0) = w(t,l) = 0,d 2 w(t,0) = d 2 w(t,l) = 0, 


w(0,x) = vvo(x),9,w(0,x) = H'l(x). 


( 1 . 1 ) 


Here w denotes the transverse deflection; d,,d x ... are the partial differentiations in t,x...\ the 
positive constants in c and D represent the unit mass, the damping coefficient and the bending 
stiffness of the panel, respectively. The axial force Q is positive or negative according to the force 
being tensile or compressive. The large panel deflection introduces an additional tension N{t) given 
by 

m « b\ o \d x w{t,x)?dx, ( 1 . 2 ) 

where b is an elastic constant. The forcing term A p denotes the pressure difference across the panal 
surfaces. The homogeneous boundary conditions mean that the panel is simply supported, and the 
initial data vv 0 and wj are given. Suppose that, without any control, the periodic solution of equation 
(1.1) excited by the pressure A p is unstable. Our problem is to stablize the panel oscillation by 
applying an appropriate control in the form of vibrational forces added to the axial force and the 
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pressure A p. 

The paper is organized as follows. To illustrate the basic ideas involved, in section 2, we 
consider the control of the Duffing equation, for which the response characteristics to a 
time-harmonic excitation is well known. The feasibility of the vibrational control can be discussed 
geometrically by referring to the response curves. Since the applicability of vibrational control is not 
limited to the structure dynamics, in section 3, a general method of vibrational control for a class of 
nonlinear evolution equations is presented. For a given unstable periodic solution, the control 
strategy is to shift the Liapunov exponent r of the vibrational equation to the negative half-line so 
that the corresponding periodic solution becomes stable. This method is applied to the nonlinear 
panel vibration problem satisfying equation (1 .1). For weak nonlinearity, analytic results are 
obtained by a perturbation analysis and the case of single-mode excitation is worked out in detail. 
Finally, in section 5, some concluding remarks are made and other possible applications such as the 
flow stabilization problem are mentioned. 

2. Control of Duffing’s Equation 

Before dealing with the nonlinear beam equation (1 .1), we consider the Duffing equation 

y + ny + Sy + Py 3 = Fcoscot, (2. 1 ) 

where the dot denotes the time derivative, the constants n,S,/3 and Fare assumed to be positive here, 
and co > 0 is the forcing frequency. For small F, by perturbation analysis [7], it is known that 
equation (2.1) has a periodic solution of the form 

y = A cos(cot + 9) (2.2) 

for some phase shift 9, where the amplitude A is related to the frequency co by the response equation 

[ 8 ], 

[(co 2 -S)A- j-pA 3 ] 2 + nco 2 A 2 = F 2 . (2.3) 

By varying the value of F, equation (2.3) yields a family of response curves in the \A\ - co plane. 
Referring to Fig.l for F = Fo,F\, the solid portion of the curve corresponds to the stable regime for 
the periodic solution (2.2), while the dotted part of the curve (between two points of vertical 
tangency) renders the solution unstable. With the frequency co fixed, point U on the Fo-curve is 
unstable, but, by changing Fo to F i, point U moves up to point S on the F’i - curve becoming a stable 
point. Therefore, in this case by adding an in-phase force with the same frequency, an unstable 
periodic motion can be stabilized. On the other hand, if the forcing amplitude F is large, the system 
may exhibit a subharmonic response. For example, consider the case of subharmonic response with 
frequency f. Again, by a perturbation analysis, it is found that equation (2.1 ) has a subharmonic 
solution of the form [8]: 


y = Acos(cot + 9\) + B cos(y/ + 0 2 ) 

where A,B and 0i,02 are the corresponding amplitudes and phases, which satisfy some response 
equations. For the subharmonics, the equation reads 

co 2 = 9a + ^y(B 2 + 2f)± [ (2^/®. ) 2 _ ^ J 2 (2 .4) 

with/= 9F/8. For F - Fq,F\ with F \ > Fo > 0, the response curves associated with (2.4) are given 
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in Fig.2. Note that for co = coq, point P on the FVcurve corresponds to a subharmonics with 
amplitude ro . However, this subharmonics will disappear when F changes from Fq to F\, since 
co < coo, and co \ is the smallest frequency for the existence of a subharmonics at F = F\. This may 
explain qualitatively why a subharmonic vibration can be suppressed in the experimental 
investigation [ 2 ] by an additive periodic force, which has the effect of changing the forcing 
amplitude F. 

In contrast with the additive vibrational control, the control can be applied parametrically. For 
instance we regard the Duffing equation (2.1) as an approximate equation for an inverted pendulum 
near the upper equilibrium position [y = 0), for which a = -8 < 0. Clearly y = 0 is an unstable 
equilibrium. If the suspension point vibrates at a high frequency v » co, the equation (2.1) should 
be replaced by [3] 


z + /id + \p(vt) - 6 ]z + fiz 3 = Fcoscot, (2.5) 

where p(z) - p{ t + 2n) is a periodic function. Without the control p, a periodic motion about y = 0 
is obviously unstable. However, by the method of averaging [3], equation (2.5) can be closely 
approximated by the averaged equation (see section 4): 

y + py + [< p 2 > -(5]y + fiy 3 = Fcoscot, (2.6) 


where< p >= 0 and 


<P n >= -57 J o P n (T)dr, for n = 1,2. 

Therefore, if a =< p 2 > -5 > 0, the periodic motion can now be stabilized as before. 

The above examples show the possibility of stabilizing periodic motions by vibrational control. 
As a generalization we consider the following control problem: 

z + pU + az + (iz 3 = j{X,cot) + h(X,vt,z) (2.7) 

where/(A,r) = r + 2n ) and h(X,r,z) = h{X, t + 2k,z) are periodic functions; A is a control 
parameter and h is a certain control function with h(0,vt,y ) = 0 and co,v are the vibration 
frequencies with v » co. The uncontrolled case corresponds to A = 0 and h = 0. Of course the 
equations (2.1) and (2.5) are special cases of (2.7). Suppose thatz = cpo(t) is unstable periodic 
solution of equation (2.7) when A = 0 and h = 0. The control objective is to choose the control 
parameter A and function h so that the corresponding periodic solution z = <p(A, /) with 
<p(0,/) = cpo(t), becomes asymptotically stable. This means analytically that the variational equation 
f° r T = (z - cp) from (2.7) has only exponentially decaying solutions. More precisely, if 

r{X, cp) = lim-f ln[y(r)|, 

then we must choose A and h such that r < 0. Obviously, unlike the optimal control, such a control, 
if possible, is far from unique. The choice of A and h, though guided by physical feasibility, is mostly 
up to the personal preference. In what follows, this control principle will be generalized to deal with 
nonlinear partial differential equations. 

3. Stabilization of Nonlinear Evolutions Equations 

In the theoretical discussion, it is convenient to consider the partial differential equations of 
interest as a nonlinear evolution equation of the form 
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(3.1) 


f ^ = B(u) + F(X,cot) + H(X,vt,u), 
\ K 0) = h. 


where u(t) is a vector in some infinite-dimensional vector space V with initial state h. The operator 
B is nonlinear, F(A, z ) = F(X,t + 2k) and H{X,t, u) = H(X,z + 2n,u) are periodic with control 
parameter A, being a scalar or vector. The control function H acts parametrically with rapid 
oscillations so that v » co. When A = 0,H(0,z,u) = 0 and the system (3.1) is uncontrolled. We are 
interested in stabilizing a periodic motion which is unstable at A = 0. If the equilibrium solution z/ 0 
of (3.1) at A = 0 is also unstable, we introduce a parametric vibrational control H to stabilize it as in 
the case of an inverted pendulum. The effect of //can be examined by the method of averaging [3], 
By a change of time from t to a = r/e with s = 1/v, the system (3.1) can be approximated by the 
averaged equation 


f ^ = B(u) + F(X,M) + D(X,u), 
1 w(0) = h. 


(3.2) 


where 


□ (A, z0 = ( H(X,vt,u )) = J Q H(X,T,u)dT. (3.3) 

The function H should be chosen so that the equilibrium solution u\ of the averaged equation (3.2) 
becomes stable. Without control (A = 0), let u = t/ /, o(0 be an unstable periodic solution of equation 
(3.1) near uq. In addition to the parametric control H, we have modulated the forcing function 
F(X,cot, ) by tuning the control parameter A so that the corresponding periodic motion satisfying the 
averaged equation (3.2) is asymptotically stable. To this end let us consider the variational equation 
of (3.2) for v = (u - y/) : 


where 


f <bL- Bx { W ,v) + Hi( X,y;,v), 
\ v(0) = g. 


(3.4) 


= B(v + if/) - B(\f/), 

//i(A,i/z,v) = D(X,v + if/)-a(X,y/), 

and g is an initial vector in V. Let ||/z|| denote the magnitude (norm) of vector h. The control 
objective is then to choose function H and parameter A in such a way that the Liapunov exponent r is 
negative, 

r(,X,H) = lim-j-ln||v(OII < 0, (3.5) 

for all g with ||g|| < 8 with some <5 > 0. For small 8, the variational equation (3.4) can be linearized 
to give 
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(3.6) 


J 

L v (°) = g> 


where A (A,/) is a linear operator defined by 

A(X,t)v = {£„[yr(A,0] + tf„[A,i/r(A,f)}v, (3.7) 

and £ u (v) = -^^-,i/ u (A,v) = — are the linearized operator of B and H at v. Note that .<4 (A, t) is 
periodic with the same period T as that of if/. Now, if K is finite-dimensional and A(X,t) is a matrix, 
then, by the Floquet theory [7], the solution of equation (3.6) can be expressed as 

v(0 = P{t)e’ R g, (3.8) 

where P{t ) = P(t + T) is a periodic matrix, and R is a constant matrix. The smallest real part of the 
eigenvalues of R yields the Liapunov exponent r. Of course the representation (3.4) holds for any 
finite-dimensional approximation of equation (3.4). Unfortunately, even the periodic function if/ is 
known, the analytical computation of the Liapunov exponent r through either (3.5) or (3.8) is 
impossible without simplifying assumptions. For example, for small amplitude vibration, the 
nonlinearity is weak so that the perturbation method and an eigenfunction expansion can be applied. 
This procedure will be illustrated in the application to the panel vibration problem. 

4. Vibrational Control of Elastic Panel 

By redefining the constants in the nonlinear beam equation (1.1) under a vibrational control, it 
yields 

d}w + pd,w - (a + P\\d x w\\ 2 )d x w + yd 4 x w = p(X,cot,x) + h(X,vt,x,w) (4.1) 

where the initial-boundary conditions are omitted, and 

IM 2 = \\dxw\ 2 dx (4.2) 

p{X,cot,x) =a p(cot,x) + (X,cot,x), (4.3) 

p i and h are the additive and parametric control forces with frequencies co » v. The physical 
constants p,P,y are positive, while a is positive or negative depending on the axial force being 
tensile or compressive. Without control, we assume that, at X = 0,p\(0,x,T ) = h(Q,o,x,w) = 0. To 
be specific, we choose the parametric control to be a vibrational axial force of the form 

h(X,vt,x,w) = q{X,vt)d 2 x w, (4.4) 

where 


q = d,q(X,vt) 

with 


Now let u\ = w and m defined by 


< q >=< q >= 0. 


(4.5) 
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ii\ = «2 + q(X, vt)d 2 x u\ 


Then equations (4.1) and (4.6) yield 


We set 


in = + ^«i) - (a + p\\d x ii\\\ 2 )dlu\ + yd x u 1 } 

- 9(9|m 2 - qd x u\) + p(X,<ot,x). 

-[!i] 


and rewrite the equations (4.6) and (4.7) in the form (3.1): 


where 


^rr = B(u) + F{X, cot) + H(X, cot, u). 


-■{pin - (a + p\\d x it]\\ 2 )dlu\ +ydju 1 } 


F{X,cot) = 


p{X,cot,*) 


H(X,vt,u) = 


z/ 1 


-pqd 2 x u\ - q 2 d x u\ - qd 2 x ii 2 


In view of equations (4.5) and (4.9)-(4. 1 1 ), by taking the time-average of equation (4.8) in o = vt 
with t = cot fixed, we get 


= B(u) + F{X,cot) + U(X,u), 


where 




< q 2 > d A x u\ 


We note that the average equation (4.12) yields a scalar equation for tv = u\ as follows, 

djw + pdjw- (a + P\\d x w\\ 2 )d 2 w+ (y +< q 2 >) = p(X,cot,x), (4.14) 

which shows that the high-frequency axial vibrational force q(X,vt ) has the effect of increasing the 
bending stiffness y by the magnitude of < q 2 > . Thus it stabilizes the system statically in general. 
Now let 
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Wo = wo(t) = 


Vo\(t) 

¥oi(t) 


be an unstable periodic solution of equation (4.8) when X = 0 (H = 0), and let 


u = i//(A,/) = 


V\(X,t) 

Vi(X,t) 


be a periodic solution of equation (4.14) with y/(0 ,t) = i/z 0 (t). Define v = (w - y/) so that v satisfies 
the variational equation (3.4). Here it can be written in the form 

=A(X,t)v + pG(X,t,v), (4.15) 

where A is a periodic linear operator and G is a nonlinear mapping defined as 


Av 


V2 

/. 



and 


f\ = -fiv 2 + (a +< q 1 > +P\\d x \f / 1 || 2 )3?vi + 2P(d x \i/i,d x v l )dly/ i - yd\v i, 
fi = ||5.vV| || 2 5ji/ / i +2(5^1,5^1)53, +115^, || 2 5 2 vi, 
with the inner product notation 

(g,h) = f g(x)h{x)dx. 

J 0 

When the nonlinear term G is dropped, equation (4.15) yields a generalized Hill’s equation, a linear 
partial differential equation with periodic coefficient: 

^f=A{X,t)v. (4.16) 

For computational purposes, introduce a complete set of orthonormal functions { e „ }, which may be 
the eigenfunctions associated with the linearized problem, or e n (x) = sm-tf-x,?} = 1,2,... By 
the expansion of the solutions of (4.12) and (4.15) into terms of e'„s as follows, 

OO ' 00 

u = ^z/„(0e„; v = Y, v„(t)e n , 

n — 1 n-\ 

their coefficients satisfy the infinite systems of coupled ordinary differential equations of the form: 

= B,{u ,«„...) + F.(X,t) + 

oo 

i= 1 > 2 , 

«-I 

where a tj = (Ae^ej). The above systems can only be solved numerically for truncated systems of 
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low dimensions. As mentioned before, if the nonlinear effect is weak, we can apply the perturbation 
analysis to approximate solutions analytically. To this end let us assume that the damping coefficient 
and the forcing amplitude are small. By proper scaling with respect to a small parameter e > 0, the 
equation (4.14) is rewritten as 

djw + Epd,w- (a + Ep\\d x w\\ 2 )d;W + ydiw = ep(X,cot,x), (4.17) 

for which we assume a > 0 and set < q 2 >= 0 for simplicity. It remains to study the problem of 
additive control. To illustrate the perturbation procedure, we will analyze the case of single-mode 
excitation in some detail. 

Let us consider the case of nth mode harmonic excitation in (4.17): 

p{X,cot,x) = F{X)sm^j-xcoscot,n = 1,2... (4.18) 

where the control parameter X modulates the forcing amplitude Fand F 0 = F( 0) is the uncontrolled 
amplitude. Then the equation (4.17) admits a single-mode solution 

w = z„(/)sin -^j-x (4.19) 


and z„ satisfies the Duffing equation: 


where 


z„ + EpQn + a n z„ + Efi„zl = eF(X) coscot, (4.20) 

[o + KTfKf-) 2 . 

A. = f (?>*■ 


The perturbation analysis of Duffing’s equation has been discussed by many authors (see e.g. [3], 
[7]). Here we adopt the method of averaging by letting 


with 


z n = 7i(0sin6J/+y/2(0 COSft, L 


(4.21) 


y i sinft)t+j> 2 cosft)t = 0, 
which are substituted into (4.20) to give 

Ti = lir[ <5 "^ 2 “ + pcoy\ + F(A)], 




(4.22) 


V. 


h = ~ +^y 2 ]. 


with 8„ = ( co 2 - a n ) and [y| 2 = y ] +y\. In the polar formal = rsin<p,y 2 = rcoscp, this equation 
becomes 


r = Ta 0 R i r >(P,(O), 
<P = 


(4.23) 


where 



R = ficor - /’simp 

_ 3 _ 
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<t> = S„- j-p n r 2 + -f cos<p. 


The solution (4.21) can be written as 


z n = r(t)cos[co/ + <p(/)]. 


(4.24) 


(4.25) 


Therefore, for z„ being periodic with frequency co,r and <p must be constants, which correspond to 
the equilibrium point (7,7p) of equation (4.23) satisfying 


j R(r,7p,co) = 0, 
\ O (7,<p,£D) = 0. 


By taking equation (4.24) into account, the above equation can be solved approximately to give 

(<5„7 - j£*7 3 ) 2 + fir 2 co 2 = F 2 (X), (4.27) 

which, by a change of notation, agrees with the response relation (2.3). Therefore for each n, the 
response curves are shown in Fig.l. Schematically, for w = 1,2,..., the response curves are plotted in 
Fig. 3. Geometrically the control strategy is to steer an unstable point U„ on the F 0 -curve to a stable 
point S„ on the F\ -curve. Analytically the stability of a periodic solution is now reduced to that of an 
equilibrium point, which can be checked more easily. To do so we form the first variational equation 
of (4.23) about ( 7,<p ) = 

f £-•£(&•/> + *,«,), 

1 , , , (4 28 ) 

where R r = d r R(r ,(p,co),R v = d v R(7 ,7p,co) and so on. Let r/(A) denote an eigenvalue ofthe 
coefficient matrix of (4.28). It can be readily verified that, by making use of (4.24) and (4.27), if 

D(X) = -£„O r ) > 0, (4.29) 

then Rep(A) < 0 so that the steady state (7,<p) is stable. This is of course the stability condition for 
the associated periodic solution. After computing the partial derivatives in (4.29), it yields 

D{X) = (j-Pn? 2 - S„) (lp n r 2 - 8 n ) + n 2 co 2 > 0, (4.30) 

with 8 n = (co 2 - a„). The above inequality determines the stability regime S in the 7 - co plane. In 
view of (4.27), 7(A) depends on the control parameter A, which will take an unstable point into the 
stable regime S. Note that from (4.30), we can get a simple sufficient stability condition: 

co 2 < a n + }P„r 2 l, (4.31) 


or 


CO 2 > a„ + \p n r 2 X. 

The above inequalities give rise to stable (shaded) sub-regions as shown in Fig.4. 
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In general all modes are excited by periodic pressure fluctuations. For instance, consider the 
harmonic forcing (4.1 8) with a general spatially dependent amplitude 

p(X,cot,x ) = F(X,x)coscot. (4.32) 

If the axial load is compressive (a < 0) and slightly exceeds the lowest buckling load y(y ) 4 , the 
parametric control q{vt) can still be used to stabilize the system statically by choosing < q 2 » |a|. 
So we remain to consider equation (4.17) by assuming a > 0 there. To apply the above perturbation 
procedure, we need to expand w in (4. 1 7) and F in (4.32) into infinite series with respect to the 
modal function e n (x ) = sin -y-x, for n = 1,2 The resulting infinite system of coupled 

nonlinear differential equations for the coefficient functions can then be treated by a perturbation 
analysis. Such a procedure developed previously for nonlinear wave equations by one of us (Chow 
[9]) can be applied here. However unlike the single mode situation, simple stability conditions such 
as (4.30) or (4.31) are no longer attainable. Though it is possible to study the stability regime 
numerically after a finite-mode approximation, this has not yet been done. 

5. Concluding Remarks. 

In the paper we present a general method of vibrational control for a certain class of nonlinear 
evolution equations with a particular reference to the nonlinear beam equation arising from the panel 
structure dynamics. The control consists of a high frequency parametric vibration and the forcing 
amplitude modulation. The high-frequency control is to affect a change in system parameter for 
static stability, while the additive control of the excitation force, if needed, is to stabilize an unstable 
periodic motion. In application to the panel structure, we show that, for a periodically excited panel 
near a buckled state, a high frequency oscillatory axial force can keep the system in the state of 
periodic motion, which can then be stabilized by an additive force modulation. The reason that we 
only control the force amplitude, instead of both the amplitude and phase is that the additive control 
is the most effective when it is in phase or out of phase with the excitation force. For a small forcing 
amplitude, a perturbation technique can be used to reduce the stabilization of a periodic motion to 
that of an equilibrium point, the latter of which is much simpler to analyze. In the case of a 
single-modal excitation, an explicit stability condition is obtained. By a finite-modal approximation, 
the stabilization problem can be studied numerically but has not yet been treated. The vibrational 
control principle described in this paper can also applied to other problems such as the flow stability 
control. Here the nonlinear evolution equation is given by the Navier-Stokes equation. For a slightly 
unstable flow, the perturbation analysis by Keller and Kogelman [10] can be employed to deal with 
the flow stabilization by vibrational control. 
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Fig. 4: Stability Regions (Shaded) for Different Modes of Excitations 
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